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ABSTRACT 

Relating the observed CO emission from giant molecular clouds (GMCs) to the underly- 
ing H 2 column density is a long-standing problem in astrophysics. While the Galactic CO-H 2 
conversion factor (Xqo) appears to be reasonably constant, observations indicate that Xco 
may be depressed in high-surface density starburst environments. Using a multi-scale ap- 
proach, we investigate the dependence of Xco on the galactic environment in numerical sim- 
ulations of disc galaxies and galaxy mergers. Xco is proportional to the GMC surface density 
divided by the integrated CO intensity, Wco, an d M^co is related to the kinetic temperature 
and velocity dispersion in the cloud. In disc galaxies (except within the central ~ kpc), the 
galactic environment is largely unimportant in setting the physical properties of GMCs pro- 
vided they are gravitationally bound. The temperatures are roughly constant at ~ 10 K due 
to the balance of CO cooling and cosmic ray heating, giving a nearly constant CO-H2 con- 
version factor in discs. In mergers, the velocity dispersion of the gas rises dramatically during 
coalescence. The gas temperature also rises as it couples well to the warm (~ 50 K) dust at 
high densities (n > 10 4 cm -3 ). The rise in velocity dispersion and temperature combine to 
offset the rise in surface density in mergers, causing Xco to drop by a factor of ~ 2 — 10 
compared to the disc simulation. This model predicts that high-resolution ALMA observa- 
tions of nearby ULIRGs should show velocity dispersions of lO^lO 2 km s — 1 , and brightness 
temperatures comparable to the dust temperatures. 

Key words: ISM: molecules - ISM: clouds - galaxies:ISM - galaxies:starburst - galaxies:star 
formation - galaxies:interactions 



1 INTRODUCTION 

Stars form in giant molecular clouds (GMCs) whose primary con- 
stituent is molecular hydrogen, H2. Because H2 lacks a perma- 
nent dipole moment, and the lowest lying excited state capable of 
quadrupole emission requires temperatures ~ 500 K to be excited, 
the physical conditions in the cold (~ 10 K) molecular gas are typ- 
ically probed via tracer molecules, rather than by direct detection 
of H 2 . 

Carbon Monoxide ( 12 C 16 0; hereafter, CO) is the second most 
abundant molecule in GMCs. Because the J=l-0 rotational tran- 
sition of CO lies only ~ 5 K above ground, has a relatively low 
effective density (~ 10 2-3 cm~ 3 ) for excitation (Evans 1999), and 
has a wavelength of ~ 3 mm which is readily observable from 
the ground, CO (J=l-0) has historically been one of the most com- 
monly used tracers of physical conditions in the molecular ISM. 



A large uncertainty in using CO to trace H2 gas is relating the 
observed CO line luminosity to the underlying H2 column density. 
However, despite the fact that CO/H 2 abundances vary strongly 
within GMCs (e.g. Sternberg & Dalgarno 1995; Lee et al. 1996; 
Hollenbach & Tielens 1999; Glover et al. 2010; Glover & Mac Low 
201 1), a multitude of observations suggests that the conversion fac- 
tor between CO and H2 is reasonably constant in Galactic GMCs, 
following the relation: 



Xco = 2 - 4 x 10 20 cm"7(K - km s" 



(1) 



where Xco is the CO-H2 conversion factor in units of H2 column 
density divided by velocity-integrated CO line intensity 1 . Lines 
of evidence for a relatively constant Xco include comparisons 
between CO luminosities and molecular column densities deter- 
mined via a variety of techniques, including dust extinction (Dick- 
man 1975), 7-ray emission (Bloemen et al. 1986; Strong & Mattox 
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Xqq is sometimes referred to in the literature as the "X-factor". We will 
use Xco an d X -factor interchangeably. 
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1996; Abdo et al. 2010b) and thermal dust emission (Dame et al. 
2001;Draineetal. 2007). 

Beyond this, the CO-H2 conversion factor appears to have the 
same relatively narrow range of values in galaxies in the Local 
Group as well (Blitz et al. 2007; Wolfire et al. 2010), though there 
may be some variations associated with metallicity (Wilson 1995; 
Arimoto et al. 1996; Boselli et al. 2002; Rosolowsky et al. 2003; 
Israel 2005; Bell et al. 2006; Leroy et al. 2006; Bell et al. 2007; Bo- 
latto et al. 2008; Leroy et al. 2011). The relatively narrow distribu- 
tion of values for Xco in Local Group GMCs 2 may arise from the 
fact that molecular clouds appear have remarkably similar physical 
properties in both the Milky Way and nearby galaxies. GMCs in the 
Galaxy and local Universe appear to have a nearly constant surface 
density of 85-100 Mq /pc 2 , obey the size-linewidth relation, and 
have relatively low kinetic temperatures of 10-20 K (Solomon et al. 
1987; Blitz et al. 2007; Heyer et al. 2009). Magnetohydrodynamic 
and radiative transfer modeling by Glover et al. (2010), Glover & 
Mac Low (2011) and Shetty et al. (2011a,b) have shown that sim- 
ulated GMCs with mean densities, sizes, velocity dispersions, and 
metallicities comparable to those found in the Galaxy naturally pro- 
duce Xco conversion factors comparable to Equation 1. 

The situation becomes more complex in starburst galaxies. By 
utilising high-spatial resolution interferometric mapping of nearby 
(ultra)luminous infrared galaxies ([U]LIRGs; L m >[10 12 ]10 11 
L ), Solomon et al. (1997), Downes & Solomon (1998) and 
Downes & Solomon (2003) have shown that the application of the 
"standard" Galactic Xco conversion factor would cause the in- 
ferred molecular gas mass to exceed the dynamical mass in these 
galaxies. In this case, the constraints on Xco are in the range 
~ 2 - 10 x 10 19 cm~ 2 /K - km/s: a factor of 2-20 lower than 
the Galactic value. Other observational evidence from local star- 
bursts (Hinz & Rieke 2006; Meier et al. 2010), the Galactic centre 
(Oka et al. 1998) and high-z submillimetre-galaxies (Tacconi et al. 
2008) have all corroborated this picture that Xco may be lower in 
regions of high molecular surface density. 

The exact origin of a lower Xco factor in starburst galaxies 
is not entirely clear. Models by Maloney & Black (1988) which 
predate the aforementioned observations, predicted that warmer 
molecular gas temperatures in infrared-luminous galaxies may 
drive a lower Xco conversion factor due to an increase in CO 
brightness temperature with kinetic temperature in optically thick 
clouds. Alternatively, Downes & Solomon (1998) suggest that the 
CO line width in starbursts traces a combination of the gaseous and 
stellar potential, rather than just the H2 mass. In the case where CO 
is optically thick, the observed velocity-integrated CO line intensity 
can increase with the velocity dispersion. Maloney & Black (1988) 
and Shetty et al. (201 lb) also postulated a similar effect if the CO 
linewidths were larger than their typical virial values. 

While the scaling of Xco with environmental parameters is 
not yet known, the ramifications are profound. For example, if Xco 
does indeed systematically vary in higher surface density environ- 
ments, our current understanding of the normalisation and index of 
the Kennicutt-Schmidt star formation rate-gas surface density rela- 
tion in star forming galaxies may change (Kennicutt 1998a; Daddi 
et al. 2010; Genzel et al. 2010). At higher redshifts, as high gas sur- 



2 -^CO a l so appears to be reasonably constant in diffuse H2 gas in the 
Galaxy (Liszt et al. 2010). This has been attributed to the offsetting effects 
of of lower CO abundances with respect to H2 (most of the carbon is in 
C+) and a large Wco/Nco ratio in low extinction gas. See the discussion 
in Pety et al. (201 1) for more details. 



face density galaxies begin to contribute substantially to the cosmic 
star formation rate density (e.g. Le Floc'h et al. 2005; Hopkins 
et al. 2010; Hopkins & Hernquist 2010), the variation of Xco with 
environment may affect observed values of the cosmic evolution 
of £^h 2 - More generally, the interpretation of forthcoming results 
from the EVLA and ALMA will be severely crippled without an 
understanding for how to relate the observed CO line flux to the 
quantity of interest: H2 gas mass. 

In this area, numerical simulations can offer some guidance. 
Indeed, understanding the origin of Xco in galaxies is difficult in 
that many of the physical parameters driving the relation are cou- 
pled. The CO-H2 conversion factor, Xco, has dimensions of: 

Xco oc E/(Wco) oc E/(a x T B ) (2) 

where Wco is the velocity-integrated CO intensity, E is the gas 
surface density, Tb is the brightness temperature of the line, and a 
is the velocity dispersion. To first order, Tb is related to the kinetic 
temperature of the gas (Tk) when the line is thermalised. However, 
the kinetic temperature of the gas can depend on molecular abun- 
dances, gas densities, dust temperatures and the background radia- 
tion field (e.g. Narayanan et al. 2006b; Krumholz et al. 201 1). The 
same physical processes which can cause changes in these param- 
eters may also cause the GMC surface densities to change as well. 
The problem is well-suited for numerical simulations. 

Building on the seminal work of Maloney & Black (1988) and 
more recent simulations of Glover et al. (2010); Glover & Mac Low 
(2011) and Shetty et al. (2011a,b), we present the first models in- 
vestigating the CO-H2 conversion factor in hydrodynamic simula- 
tions of isolated disc galaxies and disc galaxy mergers. This is the 
first paper in a series. In this work, we aim to understand whether 
Xco varies between "normal" disc galaxies and galaxy mergers, 
and if so, why. In order to do this, we couple smoofhed-particle 
hydrodynamic simulations of galaxies in evolution with dust and 
molecular line radiative transfer calculations to self-consistently 
calculate the kinetic temperature of and emissivity from GMCs in 
our models. Our main result is that higher kinetic temperatures and 
velocity dispersions in the GMCs naturally arise during mergers, 
and contribute to lower values of Xco in these systems. 

Our paper is organised as follows: In § 2, we present our radia- 
tive transfer and hydrodynamics methodology; in § 3, we discuss 
the synthetic observational and physical properties of our model 
galaxies in an effort to aid comparisons to observations; in § 4, we 
discuss how the CO-H2 conversion factor varies in disc galaxies 
and mergers in our models; in § 5, we discuss the implications of 
our findings, and in § 6 we summarise our results. 



2 METHODS 

Our goal is to simulate the emission from GMCs on galaxy-wide 
scales. This involves simulating galaxies in evolution, the physi- 
cal state of GMCs, molecular line radiative transfer through the 
clouds, and dust and molecular line radiative transfer through the 
galaxy. In this section, we describe these simulations, and the rele- 
vant assumptions that go into our modeling. This involves combin- 
ing a large number of simulation codes. In light of this, to guide the 
reader through the numerical details and equations in this section, 
we first summarise them more generally here. 

We first simulate the hydrodynamic evolution of both disc 
galaxies and mergers. It is from these simulations that we know 
the global distribution of stars, gas and metals in the galaxy, and 
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their physical properties. The radiative transfer occurs in post- 
processing. We project the physical conditions of the particles onto 
an adaptive mesh using the SPH smoothing kernel. The base mesh 
is 5 3 spanning a 200 kpc box. The cells refine recursively into 2 3 
subcells based on the refinement criteria the relative density varia- 
tions of metals {a pm / < p m >) should be less than 0.1, and the 
V-band optical depth across a cell be less than unity. The maxi- 
mum refinement level was 11, such that the smallest cells in this 
mesh are of order ~ 70 pc across. 

The surface density of and velocity dispersion within the 
GMCs are set by the physical conditions in the hydrodynamic 
galaxy evolution simulations. A subgrid prescription comes into 
play when GMCs are unresolved (i.e. when cells in the adaptive 
mesh are very large). We assume that all of the H2 mass in the 
cell is in the GMC and we calculate the HI-H2 balance via analytic 
models (described below). From this, the complete physical condi- 
tions (except for the temperature) of the GMCs are described by the 
hydrodynamic galaxy evolution simulations. The temperatures of 
the clouds are calculated by assuming thermal equilibrium between 
gas heating (by the grain photoelectric effect and cosmic rays), gas 
cooling (via molecular and atomic line cooling), dust heating (from 
the ambient radiation field), thermal dust cooling, and some energy 
exchange between gas and dust. 

With the physical properties of the galaxies and GMCs known, 
we then proceed to calculate the emergent CO emission from the 
clouds. We calculate the CO line emission from the GMCs utilising 
an escape probability formalism. The radiation from these clouds 
then interacts with other clouds in the galaxy, and the level popula- 
tions of CO are calculated by the balance of radiative absorptions, 
stimulated emission, spontaneous emission, and collisions with H2 
and He. 

At this point, the general reader should be equipped to under- 
stand the general results of this paper. For the remainder of this 
section, we elaborate on this abbreviated description. Throughout, 
we assume h = 0.7. 



2.1 Smoothed Particle Hydrodynamic Simulations of 
Galaxies in Evolution 

We simulate the hydrodynamic evolution of both idealised isolated 
disc galaxies, and mergers between these discs. The purpose of the 
hydrodynamic simulations is to calculate the spatial distribution of 
the neutral ISM, stars and metals. It is from the neutral ISM that 
we will calculate the molecular gas properties, and, as we will dis- 
cuss, the radiation from the stars and dust in the metals that de- 
termine the IR radiation field. Here, we describe the components 
of the model most pertinent to this study, namely the physics of 
the ISM and star formation prescriptions. For a more full under- 
standing of the underlying algorithms in GADGET-3, please refer to 
Springel & Hernquist (2002, 2003); Springel (2005) and Springel 
etal. (2005b) 3 . 

The galaxies are simulated with a modified version of the pub- 
licly available SPH code, GADGET-3 (Springel 2005). The ISM 
is modeled as two-phase, with cold clouds embedded in a hot, 
pressure-confining medium (McKee & Ostriker 1977; Springel & 



3 We note that Springel et al. (2005b) describes the publicly-available 
GADGET- 2, whereas the work in this paper utilises GADGET-3, a non-public 
modified version of GADGET-2. The main improvement in GADGET-3 over 
GADGET- 2 is better load balancing on parallel processors. 



Hernquist 2003). Numerically, this is realised via hybrid SPH parti- 
cles. The cold gas mass grows via radiative cooling of the hot phase, 
and cold gas is converted to hot gas through the heating associated 
with star formation. 

Stars form in the cold ISM according to a relation SFR oc 
Praid- The normalisation of this relation is set in order to match 
the local S S fr - S gas relation (Kennicutt 1998a,b; Springel 2000; 
Cox et al. 2006b). 

Supernova pressurisation of the ISM is modeled via an "effec- 
tive" equation of state (see Figure 4 of Springel et al. 2005b). Here, 
we assume a modest pressurisation of qeos = 0.25 in the Springel 
et al. (2005b) formalism. This corresponds to a mass-weighted ISM 
temperature of ~ 10 4,5 K. In the Appendix we relax the star forma- 
tion and equation of state assumptions in order to test the validity 
of our results. 

The simulations here are not cosmological: the discs are set 
up in an idealised manner in order to maximise spatial resolution. 
Here, the gravitational softening length for baryons is 100 h _1 pc , 
and 200 h~ pc for dark matter. The discs are initialised accord- 
ing to the Mo et al. (1998) formalism, and are bulgeless. They are 
embedded in dark matter halos with Hernquist (1990) density dis- 
tributions. 

In order to compare with observations in a meaningful man- 
ner, we aim to simulate galaxies comparable to those found in the 
local Universe. Accordingly, our isolated discs are initialised inside 
haloes of mass ~ 1.9 x 10 12 Mq, baryonic mass of ~ 8 x 10 10 
Mq, circular velocity of 160 km s _1 , and with 40% of the baryons 
in the form of gas. 

The mergers are binary 1:1 mergers between discs constructed 
in the same manner. We simulate three mergers of slightly higher 
mass in order to ensure that they undergo a luminous starburst 
comparable to the most extreme ones seen in the local Universe 
(~ 100 MQyr -1 ). In particular, the discs that comprise the bi- 
nary mergers have a rotation speed of 225 km s ~ , halo mass of 
~ 5 x 10 12 Mq, and baryonic mass of ~ 2.2 x 10 11 Mq. The 
mergers are set on an orbit with angles (61, </>i,#2, 4>2) = (30,60,- 
30,45), (-109,-30,71,-30) and (0,0,0,0). The angles for the first two 
orbits are arbitrary, and were chosen to represent relatively "nor- 
mal" orbits in our library of simulations. The last merger is a copla- 
nar one, and represents an extreme starburst with an extended dura- 
tion, which we include simply for comparison. We choose the first 
merger as our "fiducial" merger for the remainder of this paper as 
this particular model is well-studied in the literature 4 , and focus 
particularly on the snapshot when the star formation rate is at its 
peak. The results from all simulations are similar, and we discuss 
the minor differences that do exist when necessary. 

2.2 Physical Properties of Giant Molecular Clouds 

We assume that the entire neutral mass in a given cell is locked in 
a cloud which is spherical, isothermal, and of constant density. We 
determine the surface density of the neutral gas via 

Edoud = max(E ce ii, 100 M pc" 2 ) (3) 

where E co n is the surface density of the cell in the SPH simulation. 
In this model, when the cloud is resolved, we use the surface den- 
sity as calculated in the simulations. When the cloud is unresolved, 
we adopt a subresolution surface density comparable to observed 

4 In § 3, we discuss the physical and simulated observational properties of 
our fiducial merger to highlight its similarity to observed local galaxies. 
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values of GMCs (e.g. Solomon et al. 1987; Blitz & Rosolowsky 
2006). 

We then determine the H2 fraction of the neutral ISM utilising 
the analytic formalism of Krumholz et al. (2008, 2009a) and Mc- 
Kee & Krumholz (2010). This prescription aims to model the bal- 
ance between the dissociation of molecules by Lyman- Werner band 
photons, and the formation of molecules on dust grains. We refer 
the readers to the aforementioned papers for the full derivation, and 
simply repeat the numerical prescription here. The molecular frac- 
tion is given by: 



1-?- 



' 4 1 + 0.25s 
for s < 2 and /m = for s > 2. s 



(4) 



O.Olx )/(0.6t c ), where x 



= ln(l + 0.6 X + 
0.76(1 + 3.1Z' ' 365 ), and r c = 



O.O66E c ioud/(M pc~ 2 ) x Z'. Z' is the metallicity divided by the 
solar metallicity. This formalism for deriving /h2 assumes chemi- 
cal equilibrium. 

It is worth a quick note that there are numerous prescriptions 
for determining the HI/H2 balance in the ISM of simulations, some 
of which include time-dependent chemistry. Blitz & Rosolowsky 
(2006) developed an empirical pressure-based methodology for 
calculating the H2 fraction in the neutral ISM, based on obser- 
vations of local galaxies. Similarly, both semi-analytic models 
(Obreschkow et al. 2009; Obreschkow & Rawlings 2009), as well 
as full numerical solutions exist which model the effect of disso- 
ciating photons through models of galaxies (e.g. Pelupessy et al. 
2006; Dobbs et al. 2008; Robertson & Kravtsov 2008; Pelupessy & 
Papadopoulos 2009; Gnedin et al. 2009; Gnedin & Kravtsov 2010). 
We motivate our usage of the analytic prescription of Krumholz 
et al. (2009a) for two reasons. First, some observational evidence 
suggests that on small scales (< 100 pc), Equation 4 may fare 
better than pressure-based prescriptions in describing the state of 
the neutral ISM in low-metallicity dwarf galaxies (Fumagalli et al. 

2010) . Second, a comparison between Equation 4 and a numerical 
treatment of time-dependent chemical reaction network and radia- 
tive transfer in galaxies suggests that the analytic approximation 
is reasonable at metallicities above 0.01 Zq (Krumholz & Gnedin 

2011) . Because we aim to model actively star-forming systems in 
this work, we find that the mass-weighted metallicity of our model 
clouds is always higher than this fiducial value and expect that the 
analytic approximation is therefore reasonable. 

With E c ioud and Mh2 defined, the radius of the cloud is 
known. In order to account for the turbulent compression of gas, 

2 /2 

we scale the volumetric densities of the GMCs by a factor e />' 
where numerical simulations show 

a 2 » ln(l + 3M! D /4) (5) 

where Mid is the 1 dimensional Mach number 5 of the turbulence 
(Ostriker et al. 2001; Padoan & Nordlund 2002, see also Lemaster 
& Stone (2008)). Because the temperature calculation is dependent 
on the density of the GMC (see below), solving for the density and 
temperature simultaneously is a computationally lengthy process 
for the multi-million-cell grids that concern us. Thus, to calculate 
the turbulence-driven density enhancement, we assume the temper- 
ature of the GMC is 10 K, which as we shall show, is a good ap- 
proximation for the bulk of the GMCs in these simulations. 



We calculate the ID velocity dispersion in the cloud: 

a = max(i7 C eii,0Vir) 



(6) 



where <7 C eii is the mean square sum of the subgrid turbulent ve- 
locity dispersion within the GMC and the resolved nonthermal ve- 
locity dispersion. The subgrid turbulent velocity dispersion is cal- 
culated from the external pressure from the hot ISM (Robertson 
et al. 2004) using a 2 — Pj p ce u though we impose a ceiling of 10 
km s _1 which comes from average values found in turbulent feed- 
back simulations (e.g Dib et al. 2006; Joung et al. 2009; Ostriker 
& Shetty 201 1). The resolved nonthermal component is calculated 
by finding the turbulent velocity dispersion of the nearest neigh- 
bouring cells in the simulation. In detail, we calculate the standard 
deviation of the velocities of the nearest neighbour cells in the x, y 
and z directions, and define the nonthermal velocity dispersion as 
the mean of these. In cases where the GMC is unresolved, a floor 
ovir is set by assuming the cloud is in virial balance with a virial 
parameter a v ir = 1, for a v - lr = 5a 2 ir R/(GM), so that 



2.2 km s" 



M 



10 5 M 



1/4 



(7) 



for E c i ud = 100 Mq pc -2 where M is the mass of the cloud. 

Finally, we calculate the temperature of the model GMCs. The 
model is based on that developed by Krumholz et al. (2011), and 
we describe the relevant details here as it is an important aspect of 
our model. The temperature of the molecular ISM is determined by 
a balance of heating and cooling processes in the gas, heating and 
cooling of the dust, and a dust-gas thermal exchange. For the gas, 
we consider grain photoelectric heating at a rate per H nucleus r pe , 
cosmic ray heating at a rate Tor, and cooling via either CII or CO 
line cooling at a rate Ai ine . The dust can be heated by the back- 
ground infrared radiation field at a rate Fdust, and cool via thermal 
emission at a rate Ad us t- Finally, there is an energy exchange be- 
tween dust and gas at a rate vfgd where V&gd is positive if the dust is 
hotter than the gas. If the gas and dust are in thermal balance, then 
we have the following equations: 



Tpc + TCR - Ali ne + *gd = 

Tdust — Adust — * g d = 



(8) 
(9) 



The equation is solved by simultaneously iterating on the tempera- 
tures of the gas and dust 6 . 

The grain photoelectric heating rate is assumed to be attenu- 
ated by half the mean extinction of the cloud (as the heating rate is 
expected to decrease toward the cloud interiors) and is given by: 



erg s 



(10) 



where G' is the FUV intensity relative to the Solar neighborhood, 
and <Td is the dust cross section per H atom to UV photons. Here, 
we assume that the Go = 1 and Od = 1 x 10~ 21 cm~ 2 . Test 
models in which we scale Go by the star formation rate density 
compared to that found in the solar neighbourhood (e.g. Ostriker 
et al. 2010) have similar results to those presented in this work, and 
are presented in the Appendix. 

The cosmic ray heating rate is given by: 



Tor = C'<?cr s 1 



(11) 



5 We note that other authors have found a range of possible forms for 
Equation 5. For example, Lemaster & Stone (2008) find cr 2 « 0.61n(l + 
0.5Mf D ), while Price et al. (2011) find cr 2 « ln(l + i/9(M 2 D )) where 
M3D is the 3D Mach number. 



6 We note that this dust temperature is not always the same as the tempera- 
ture calculated by SUNRISE (§ 2.3). This makes little difference on the final 
results. We discuss this in more detail in the Appendix. 
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where £' is the cosmic ray ionisation rate (here, assumed to be 2 
x 10~ 17 Z's^ 1 ), and gcR is the thermal energy increase per cos- 
mic ray ionisation. For H2, qcn ~ 12.25 eV (though note that 
this value is quite uncertain; see discussion in Appendix A4 of 
Krumholz et al. 2011), and for HI, gci) = 6.5 eV (Dalgarno & 
McCray 1972). We utilise a constant cosmic ray heating rate for all 
simulations. Some models suggest that there may be enhanced cos- 
mic ray fluxes during starbursts which would increase the H2 gas 
temperature (Papadopoulos 2010; Papadopoulos et al. 2010), and 
further enhance the effects found in our Results section. 

Finally, in a subset of models we have explored the poten- 
tial effects of turbulent heating on molecular clouds. In unresolved 
GMCs we can estimate this heating rate based on numerical exper- 
iments on the rate of turbulent dissipation: r tur b ~ 1-5 x a 3 /R, 
where R is the GMC radius (McKee & Ostriker 2007). For re- 
solved GMCs, we can measure the turbulent heating rate directly 
from the code. Bulk turbulent motions can be converted to heat 
through two pathways: adiabatic compression and viscous dissi- 
pation. The compressive heating rate per unit mass is r com p = 
P(V • v)/p, and we can evaluate this directly from the density 
and velocity fields output by Gadget. The viscous dissipation rate 
per unit mass is r v i sc = (7r v i sc • V) • v/p, where 7r v i sc is the vis- 
cous stress tensor. The code relies on implicit dissipation rather 
than an explicit viscosity, but we can estimate the viscous heat- 
ing rate produced by that implicit dissipation by noting that the 
Reynolds number must be ~ 1 on the resolution scale of the 
code (Offner et al. 2009). This implies that the dynamic viscos- 
ity is r) ~ pvh, where h is the SPH smoothing scale. Given this 
approximation, the components of the viscous stress tensor are 
TiYj.visc = i][dvi/dxj +dvj/dxi — (2/3)dvi/dxjSij], and which 
we can again evaluate directly from the density and velocity fields 
output by Gadget. We find that the effects of turbulent heating are 
modest in both the resolved and unresolved cases. In our fiducial 
merger including viscous dissipation reduces Xco by ~ 30%, 
while the the fiducial disc it reduces Xco by less than a few per- 
cent. Hereafter we neglect this heating term, though we note that 
including it would only enhance the results we present below. 

The line cooling is assumed to occur via either CII or CO 
emission. The fraction of hydrogen for which the carbon is mostly 
in the form of CO is well-approximated by the following re- 
sult from both semi-analytic (Wolfire et al. 2010) and numerical 
(Glover & Mac Low 201 1) work: 

G' 

-4(0. 53-0. 0451n 2-_7r -0.0971nZ')/A v 

fco = /H2 X e "H/cm J (12) 

When this fraction is above 50%, we assume the cooling hap- 
pens predominantly via CO line cooling; else, the cooling occurs 
via CII emission. The cooling rate is calculated via an escape prob- 
ability formalism utilising the public code of Krumholz & Thomp- 
son (2007). We describe the equations for the line radiative transfer 
(both within clouds, as is pertinent to calculating the cooling rates, 
and across the model galaxy, in § 2.4). 

The dust cooling rate is: 

A dust = K{T d ) m caTl (13) 

We assume the bulk of the dust heating happens via IR radiation 
as IR radiation likely dominates the heating over UV flux in the 
optically thick centres of GMCs. The IR radiation field is known 
from SUNRISE dust radiative transfer calculations (which will be 
described in § 2.3). 

Finally, the dust and gas exchange energy via: 

* gd = a gd n H T g 1/2 (T d -T g ) (14) 



where the thermal gas-dust exchange rate is a gd = 3.2 x 10 34 Z' 
erg cm 3 KT 3/2 for H 2 , and a gd = 1 x 10~ 33 Z' erg cm 3 KT 3/2 
for HI (Goldsmith 2001). 

2.3 Dust Radiative Transfer 

In order to calculate the background radiation field from stars and 
the dust temperature, we perform dust radiative transfer calcula- 
tions with the publicly available code SUNRISE. A full description 
of the algorithms can be found in Jonsson (2006); Jonsson et al. 
(2010), and Jonsson & Primack (2010). Here, we summarise the 
aspects of the simulations most relevant to this study. 

The sources of radiation in the model galaxies are stellar clus- 
ters and accreting black holes. The stellar clusters emit a template 
spectrum derived from STARBURST99 calculations, with the metal- 
licities, masses and ages known from the GADGET-3 simulations. 
The AGN emits a spectrum based on observations of unreddened 
quasars (Hopkins et al. 2007b), though has little effect in the calcu- 
lations here (see the Appendix). 

The substructure of the ISM on scales below the smoothing 
length of the SPH simulations is unresolved. We assume that star 
clusters with ages < 10 Myr reside in natal birthclouds, and mod- 
ulate their SED accordingly. These birthclouds contain HII regions 
and photodissociation regions (PDRs) whose SEDs are calculated 
utilising ID MAPPINGSIII photoionisation models (Groves et al. 
2004, 2008; Jonsson et al. 2010). The time-averaged PDR cover- 
ing fraction is a free-parameter. We assume a constant fraction of 
/pdr = 0.3, corresponding to a covering lifetime of ~ 2 — 3 Myr 
(Groves et al. 2008). This value is motivated in part by simulations 
by Jonsson et al. (2010) which showed covering fractions compara- 
ble to these result in synthetic SEDs of disc galaxies comparable to 
the SINGS sample (Kennicutt et al. 2003). Changing this parameter 
has minimal effects on the final results of this paper: we quantify 
this and other potential effects of the subresolution modeling in the 
Appendix. 

When radiation leaves either the naked stellar cluster (with age 
> 10 Myr), or the HII region/PDR (for younger clusters), it is al- 
lowed to interact with the diffuse ISM. We assume the remaining 
cold molecular phase has an negligible cross-section for interac- 
tion, though test the effects of this assumption in the Appendix. 

The dust mass in the diffuse ISM is calculated assuming a con- 
stant dust to metals ratio of 0.4 (Dwek 1998; Vladilo 1998; Calura 
et al. 2008), where the metallicity distribution is known from the 
SPH calculations. We use the Weingartner & Draine (2001) dust 
model with R = Av/Eb-v = 3.15, as updated by Draine & Li 
(2007). The dust and radiation field are assumed to be in radiative 
equilibrium, utilising the methodology of Juvela (2005) for calcu- 
lating the converged radiation field. When the radiation field has 
converged, we calculate the dust temperature in each cell by iterat- 
ing equations 6-8 of Jonsson & Primack (2010) utilising a Newton- 
Raphson scheme. 

2.4 Molecular Line Radiative Transfer 

Finally, with information about the spatial distribution of GMCs in 
the model galaxies, and their mean H2 fractions, densities, temper- 
atures, velocity dispersions and kinematics through the galaxy, we 
are prepared to calculate the emergent CO line emission from the 
model galaxy. This involves two stages. First, we calculate the es- 
cape probabilities of the CO lines from the GMCs. We then track 
the propagation of these photons through the model galaxy as they 
potentially interact with other GMCs. 
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Generally, CO line emission is set by the level populations. 
The source function from a given region for a given transition from 
upper level to lower level u — > I is given by: 



S v = 



(15) 



(niBiu - n u B ut ) 

where A u i, Bi u and B u i are the Einstein coefficients for sponta- 
neous emission, absorption, and stimulated emission, respectively, 
and n are the absolute level populations. 

We first calculate the level populations within and probability 
for photons to escape from the individual GMCs in the galaxy util- 
ising the publicly available code described in Krumholz & Thomp- 
son (2007). The levels are assumed to be in statistical equilibrium 
and determined through the rate equations: 



5^(Cj„ + A„A,„)/i = 



^2(C ut + /3 ui A ul ) 



fu 



(16) 
(17) 



where C are the collisional rates, / the fractional level populations, 
and /3ui is the escape probability for transition u — ¥ I. The rate 
equations can be rearranged as an eigenvalue problem, and solved 
accordingly. 

The escape probability, /3 u i can be approximated by by relat- 
ing it to the optical depth in the line, t„j (Krumholz & Thompson 
2007): 



(18) 



l + 0.5r ui 

In the escape probability formalism, the optical depth of the 
line through the cloud can be represented as: 



Tul 



gi 16(27r) 3 / 2 cr 



Urn 



(19) 



where Q is the abundance of CO with respect to H2, gi and g u are 
the statistical weights of the levels, Nh2 is the column density of 
H2 through the cloud, X u i is the wavelength of the transition, and 
a is the velocity dispersion in the cloud. Equations 16-19 are it- 
erated upon utilising the Newton-Raphson Method until the escape 
probabilities and level populations within the GMCs are known. 

With p u i calculated, we determine the effects of radiation 
from individual GMCs on other GMCs in determining the final 
level populations utilising the 3D non-LTE Monte Carlo radiative 
transfer code TURTLEBEACH (Narayanan et al. 2006b, 2008b). We 
begin with the level populations found from the escape probabil- 
ity calculations as a guess, and emit model photons from each 
GMC isotropically with direction drawn randomly, and emission 
frequency drawn from a Gaussian profile function: 



v 



v ■ 11— ) /Avd 



(20) 



where vq is the rest frequency of the line, v is the velocity of the 
cloud in the direction of the photon's emission, c is the speed of 
light, and A^d is the doppler-width of the emission line. 

When the photon passes through a cell, it interacts with a 
GMC and sees an opacity of: 



a"' (gas) = Vfln 



47T 



■cj>{v){mBi u - riuBui) 



(21) 



where Vmi is the volume filling factor of the spherical GMC. We 
neglect absorption by dust in this model. 



After all GMCs have emitted some number of model photons, 
the level populations in the GMCs are updated by assuming de- 
tailed balance: 



ni 



^2l3ikAi k + ^2(B lk J„ +Ci k ) = 

k<l kjtl 

^^n k fikiAki + ^^n k {B k iJ v + Cm) 



(22) 



k>i 



k^i 



where Cik and Cki are the collisional rates, and f) only exists for 
transition k — > I such that k = I + 1. Equations 22 are solved via 
Gauss-Jordan matrix inversion. 

This process is iterated upon until the level populations have 
achieved convergence. Here, we demand that they not vary by more 
than a fractional difference of 1 x 10~ 3 for at least 3 iterations. 

Once the level populations have been solved for, we build the 
formal spectrum by choosing an (arbitrary) viewing angle, and in- 
tegrating along lines of sight (e.g. Walker et al. 1994): 



,M1 



-(tot) 



(23) 



Tests of TURTLEBEACH against the publicly available Leiden 
Benchmarks (van Zadelhoff et al. 2002) are presented in Narayanan 
et al. (2006b). We obtained our coefficients from the Leiden Atomic 
and Molecular Database (Schoier et al. 2005). We assume a frac- 
tional carbon abundance of 1.5 x 10~ 4 , though the abundance of 
CO with respect to H2 is given by Equation 12. 



3 OBSERVATIONAL AND PHYSICAL PROPERTIES OF 
SIMULATED GALAXIES 

As we aim to compare potential variations in Xco in our simu- 
lated galaxy mergers to those that are actually observed, it is worth 
briefly comparing the physical and synthetic observational prop- 
erties of our model galaxies to real galaxies. Our fiducial merger 
has been well-studied in the literature, and is very much an aver- 
age merger simulation as far as the range of simulated SFRs, black 
hole accretion rates and bolometric luminosities. While the pro- 
cesses described in this section generically describe gas-rich merg- 
ers, what we summarise here has been calculated and published 
previously explicitly for our fiducial model. 

The merger goes through elevated star formation rate upon 
first passage as tidal torques on the gas cause the gas to lose an- 
gular momentum and fall toward the centres causing high-density 
regions (Mihos & Hernquist 1996; Hopkins et al. 2006; Narayanan 
et al. 2008c; Juneau et al. 2009). The galaxy undergoes a starburst 
upon final coalescence, and peaks in its bolometric luminosity. Ra- 
diative transfer post-processing on these models have found that the 
model galaxy is then visible as a ULIRG (Chakrabarti et al. 2007). 
The same gaseous inflows can drive sufficient black hole growth 
to result in optical quasar activity (Di Matteo et al. 2005; Hopkins 
et al. 2005b,a, 2006). Effects of the AGN feedback can be seen in 
both the warm infrared colours of the galaxy (Younger et al. 2009), 
as well as molecular outflows (Narayanan et al. 2006a, 2008b). 

The truncation of the starburst by a combination of gas con- 
sumption and AGN feedback can render the galaxy observable as 
an E+A post-starburst (Snyder et al. 2011) before it evolves into 
into a dead early-type (Springel et al. 2005a; Hopkins et al. 2007a) 
with colours comparable to those observed on the red-sequence 
(Springel et al. 2005a; Hopkins et al. 2008b,c). The kinematic (Cox 
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Figure 1. Velocity integrated brightness temperature Wco maps for model disc galaxy (left), and fiducial merger snapshot (right). The colour scales correspond 
to the colour bars on the right of each panel, and the units are logio (K-km s — 1 ). The centre panel shows the model disc again, but with the same colour scale 
as the merger for comparison. Low gas kinetic temperature (~ 10 K) and velocity dispersions cause the bulk of the disc galaxy to have intensities of ~ 10 
K-km s -1 . In contrast, the velocity dispersion within GMCs in the ULIRG can be many tens of km s _1 , with gas kinetic temperatures near 50 K. Summed 
over a sightline, the observed gas intensity can be > 10 3 K-km s — 1 . 



et al. 2006c), X-ray (Cox et al. 2006a), nuclear emission (Hopkins 
et al. 2008d, 2009), and molecular disc properties (Xu et al. 2010) 
of this merger remnant have all been studied and found to be com- 
parable to those observed. Similarly, the remnant lies on the funda- 
mental plane (Di Matteo et al. 2005; Robertson et al. 2006; Hopkins 
et al. 2008a). 



4 RESULTS 

4. 1 GMCs in "Normal" Discs 

In the far left panel of Figure 1, we show the velocity-integrated 
brightness temperature map of the model disc galaxy. As expected, 
the central regions are the brightest, and the outer disc has little CO 
emission. In the top left panel of Figure 2, we plot the emission- 
weighted distribution of Xco values for the GMCs in our model 
disc galaxy and the fiducial model merger. We additionally plot 
the distribution of GMC physical properties in both the disc and 
merger. We will return to this plot frequently throughout this sec- 
tion and the next. 

The luminosity- weighted Xco in °ur model disc is ~ 4 x 10 20 
cm _2 /K-km s -1 with a relatively narrow dispersion. The disper- 
sion is narrow because the surface densities, kinetic temperatures 
and velocity dispersions of the model disc GMCs show fairly little 
variation. To remind the reader, the column densities in the GMCs 
in our disc galaxy are set to be the surface density of cold gas in 
the cell. When the GMC is unresolved in the simulation, we set the 
subgrid value of the surface density to E c i ou d = 100 Mq /pc 2 . 
This value was chosen to match the roughly constant surface den- 
sity of Galactic molecular clouds. Nearly all of the GMCs in the 
model disc take on this value for a surface density. 

The kinetic temperatures of GMCs in the disc have a relatively 
tight distribution near 10 K, as shown in Figure 2. Because the 
GMCs have a relatively low density compared with starbursts (the 
mass-weighted value is ~ 500 cm -3 ), there is little coupling with 
the dust grains (which are a factor of a few hotter; Figure 2). Thus 
the temperature is primarily determined by molecular line cooling, 
and heating by cosmic rays and the grain photoelectric effect. The 
kinetic temperature helps to set the brightness temperature, though 
the two are not identical. The emission-weighted brightness tem- 
perature for the merger(disc) are ~ 50(7)K. 



Finally, the distribution of velocity dispersions in the GMCs 
is fairly narrow. Recalling § 2, the velocity dispersion of the clouds 
is taken by calculating the dispersion amongst the cell's nearest 
neighbours, with a subgrid model for unresolved clouds (Equa- 
tion 7). Because the disc is dynamically cold, the velocity disper- 
sions are primarily set by the latter case. This results in an emission- 
weighted velocity dispersion within GMCs in the model disc of ~ 3 
km s _1 , with a maximum of ~ 15 km s . These values compare 
favourably with the velocity dispersions reported in the comprehen- 
sive survey of Solomon et al. (1987), and the more recent review 
by Blitz et al. (2007). 

We can ask why the simulated Xco from the model galaxy 
is comparable to the Galactic average, Xco ~ 2 — 4 x 
10 20 cm -2 /K-km s _1 . In principle this occurs because the physi- 
cal conditions in the model GMCs by and large match those of 
observed GMCs in the Milky Way. In this sense, the fact that our 
model value for Xco in quiescent discs matches that of the Galaxy 
is by construction. However, there are two salient points here. 

First, it is important to remember that we allow for the possi- 
bility that the galactic environment can set the physical conditions 
in the GMCs if the pressure is sufficiently high. The fact that the 
default value for the surface density and velocity dispersions in the 
clouds is typically used is a statement that the galactic environment 
in the model disc galaxy is not sufficiently extreme to cause sig- 
nificant changes in the surface densities, temperatures, or velocity 
dispersions in the GMCs from the Galactic values. As we will see 
in the subsequent section, this is not the case in mergers. 

Second, the subresolution values for the GMCs are not with- 
out physics. GMCs through the Local Group are observed to obey 
the Larson (1981) relations: they follow a linewidth-size relation- 
ship with a oc R 5 , they have virial parameters a ~ 1, and they all 
have roughly the same surface density £ ~ 100 Mq pc -2 (Blitz 
et al. 2007; Bolatto et al. 2008). The origin of these observed rela- 
tionships is debated, but their universality argues for some sort of 
internal regulation mechanism operating in GMCs (e.g. Krumholz 
et al. 2006; Shetty & Ostriker 2008). Regardless of the underlying 
mechanism, though, our subgrid model is not simply tuned to re- 
produce the "right" Xco - Instead, it models the real physical prop- 
erties of GMCs 7 . 



We note, however, that clouds need not be virialised to have Xco com- 
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Finally, it is important to note that the results presented in this 
section do not necessarily translate to disc galaxies at high-redshift. 
Galaxies on the "main- sequence" of star formation rates at high-z 
(e.g. Noeske et al. 2007a,b) still form stars at rates comparable to 
present-epoch mergers (Daddi et al. 2005, 2007) because they have 
very high gaseous surface densities, though they may be morpho- 
logically classified as discs (Forster Schreiber et al. 2009). 

4.2 Xco in Merger-Driven Starbursts 

We now turn to Xco in galaxy mergers. Before embarking on the 
remainder of this section, it is important to emphasise that the many 
of the GMCs in the model starburst are resolved (this is clear from 
the Sh2 panel in Figure 2). Thus, the derived values for Xco are 
independent of subresolution assumptions. 

During the merger, gas is funneled toward the nuclear regions, 
causing dense concentrations of molecular gas (Barnes & Hern- 
quist 1991, 1996). The surface densities of the GMCs in our simu- 
lations rise accordingly. In principle, this would cause a rise in the 
CO-H2 conversion factor (c.f. Equation 2). However, during the 
merger-induced starburst, the increase in velocity-integrated line 
intensity exceeds the rise in surface density, causing Xco to drop 
from the Galactic value. 

In Figure 3, we show the evolution of the star formation rate, 
gas temperature, velocity dispersion and Xco as a function of time 
for the three model galaxy mergers. The shaded region denotes the 
range of mean values among the GMCs within the merger models 
at each timestep, i.e. at time the lowest point outlined in gray cor- 
responds to the lowest galaxy-averaged value of the three merger 
models, and the highest point in gray corresponds to the highest 
galaxy-averaged value among the three. The time axes are centred 
around the peak in the starburst for each model. When the galaxies 
merge, the discs are destroyed. During this time, the dominant con- 
tributor to a within the GMCs is the nonthermal component derived 
from the local resolved velocity dispersion of the gas. The nonther- 
mal velocity dispersion is driven by the dynamics during the galaxy 
merger and mixing of stellar mass with the H2 gas. 

During final coalescence in the merger, when the SFR peaks 
at a few hundred MQyr -1 , the fraction of dense gas rises, a result 
verified both in theoretical models (e.g. Mihos & Hernquist 1994, 
1996; Narayanan et al. 2008a,c; Bournaud et al. 2011), and obser- 
vations (Juneau et al. 2009). The mass-weighted mean GMC den- 
sity rises to roughly > 10 4 cm~ 3 , compared to ~ 500 cm -3 in the 
model disc. At these high densities, the energy exchange between 
dust and gas becomes efficient, and the gas temperatures begin to 
approach the dust temperatures. At the same time, the dust is being 
heated by an amplified radiation field due to the merger-induced 
starburst. This is demonstrated explicitly in Figure 2, where we 
show the dust and gas (kinetic) temperature distributions of the 
GMCs in the model galaxies. The mean gas temperature is higher 
by a factor of a few than the roughly ~ 10 K GMCs in the model 
disc. The rise in the gas kinetic temperature during the starburst is 
shown in Figure 3 as well 8 . 

parable to observed galactic values. Provided that E c [ ou( j, T, and a remain 
within a modest range of values, Xco ~ 2 - 4 X 10 20 cm~ 2 /K-km s — 1 
(Shettyet al. 2011b). 

8 We remind the reader that we adopt a constant cosmic ray ionisation 
rate in all models. If cosmic ray energy densities increase in starburst en- 
vironments as suggested by recent observational (Abdo et al. 2010a) and 
theoretical (Papadopoulos 2010; Papadopoulos et al. 2010) work, then the 
gas temperatures would further increase, causing Xco in starbursts to fall 



The large molecular gas densities in the merger also mean the 
CO is thermalised in the ground state transition. When level popu- 
lations are in local thermodynamic equilibrium (LTE), their source 
function can be described by the Planck Function. In this limit, 
the source function rises with temperature. Hence, the rise in gas 
kinetic temperature during the burst contributes to driving Xco 
down. 

The combination of the increased velocity dispersion and the 
brightness temperature combine to exceed the increase in surface 
density, which causes a depressed mean Xco from the Galactic 
value during the merger. We now return to the first panel Figure 2 to 
explicitly compare Xco in the merger against the disc galaxy. We 
see that Xco has a broad distribution for the model merger. A num- 
ber of GMCs outside of the nucleus are similar to the disc galaxy's 
in terms of their physical properties. These GMCs are unresolved in 
our simulations (owing to the fact that they reside in lower-density 
environments), and thus take on surface density and velocity disper- 
sion values comparable to those observed in the Galaxy. However, 
the GMCs toward the centre of the galaxy all have larger surface 
densities, velocity dispersions, and kinetic temperatures. The lat- 
ter two combine to depress Xco compared to the values seen in 
the disc by a factor of ~ 5 — 10. Because most of the mass in 
the merger is in the central regions, the luminosity-weighted mean 
is low. Test simulations with fixed temperatures or velocity disper- 
sions show that the increased kinetic temperature and velocity dis- 
persion in the gas contribute roughly equally to the increased line 
intensity in the merger simulation. This is somewhat apparent from 
Figure 2, where we see similar distributions values for the kinetic 
temperatures and velocity dispersions in the gas. 

The magnitude by which Xco decreases is dependent on the 
strength of the merger. Turning to Figure 3, we see a range in Xco 
values during the burst. The model with the largest Xco during 
the burst corresponds to the lowest peak SFR. A key point of this 
aspect of the model is that there is no "merger value" of Xco : Xco 
depends on the physical parameters of the emitting galaxy. 

What happens in the post-starburst stage is also highly merger- 
specific. During this phase, the galaxy is a gas-poor early type. Gen- 
erally, the gas has a large velocity dispersion for at least a dynam- 
ical time after the burst. This is consistent with what was seen in 
simulations of CO gas in high- z submillimetre galaxies (Narayanan 
et al. 2009). During this phase, it is less trivial to simply relate the 
observed Xco to the gas velocity dispersion and temperature as 
there is a much larger dispersion in molecular gas fractions and 
CO abundances. This owes to the fact that there are highly vary- 
ing physical conditions in the post-burst galaxy, which drive strong 
variations in the H2 and CO abundances. Some of the simulations 
return to a Galactic Xco value quickly, while others remain low. 

In summary, during the merger-induced starburst, Xco drops 
in galaxies from the standard Galactic value due to increased gas 
temperatures and velocity dispersions. During this time, the CO 
abundances are ~ 1 x 10~ 4 /H2 and molecular gas fractions near 
unity in the main CO emitting region. In the resulting gas-poor 
merger remnant, the dynamical and thermal history can vary from 
model to model, and the evolution of Xco is less uniform among 
mergers. 



even further (though see Bell et al. 2006, for an expanded study on the role 
of cosmic rays in the X-factor at various extinctions). 
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Figure 2. Emission-weighted distributions of Xco values and physical properties for the GMCs in the model disc galaxy (red-dashed line) and fiducial 
merger model (solid-black line). The ordinate values are normalised. Starting from the top left, and going clock-wise, the plots show Xqq, velocity dispersion, 
kinetic temperature, brightness temperature, dust temperature, and GMC surface density. The disc galaxy predominantly has GMCs with physical properties 
comparable to the Milky Way's, and thus has similar Xqq values. The GMCs in the merger show a broad distribution in Xco values, with a lower mean than 
the disc galaxy. The lower Xco owes to larger gas temperatures (which are larger due to efficient coupling with the warm dust at high densities) and large 
velocity dispersions in the gas. The numbers in each panel refer to the logio of the emission-weighted mean value, and the black (top) number corresponds to 
the merger whilst the red (bottom) corresponds to the disc. Because the numbers correspond to the logio of the mean in the physical quantities, they will have 
larger values than one would pick by eye in the log-log plots. 



4.3 The Variation of Xco with Galactocentric Radius 

With the concepts presented in § 4.2, we are now in a position to 
understand how Xco varies in galaxies as a function of spatial lo- 
cation. In Figure 4, we show the simulated Xco maps for the disc 
galaxy and fiducial merger, and in Figure 5, we plot the values for 
Xco in the GMCs in our model disc galaxy and fiducial merger as 
a function of radius from the centre of the galaxy. The Xco values 
from Figure 5 come from the map in Figure 4. The Xco values 
are binned in bins of distance, and represent the emission-weighted 
mean within a given distance bin. The bars denote the range of 
Xco values seen in a given distance bin. 

Xco in the centre of the model disc galaxy is systematically 
lower than in the rest of the galaxy. In particular, a number of GMCs 
along the line of sight have velocity dispersions larger than the typ- 



ical virialised values, with values elevated by a factor of ~ 2. Sim- 
ilarly, due to the elevated densities in the nucleus combined with a 
warmer dust temperature, the gas temperatures of some GMCs can 
reach values up to 15 K. This causes Xco in the central kiloparsec 
to generally display the lowest values in the galaxy. Depressed val- 
ues of Xco from the Galactic mean have been observed in at least 
a few GMCs toward the Galactic Centre (Oka et al. 1998). It is im- 
portant to note that the regions where Xco > 10 21 cm -2 /K-km s -1 
represents much of the area, but a negligible fraction of the gas 
mass in the galaxy. This is evident from Figure 2. 

In the fiducial model merger, unlike the situation with the 
model disc galaxy, we see no clear trend in Xco with galacto- 
centric radius. Due to the violent nature of the gasdynamics dur- 
ing the merger, gas of a variety of physical conditions is mixed 
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Figure 3. The evolution of the star formation rate, emission-weighted velocity dispersion, kinetic temperature, and Xqq of the individual GMCs within the 
galaxies for all three merger models as a function of time. The time axes are centred for each model about the point of maximum star formation rate. The grey 
shaded region denotes the range in emission-weighted mean values for all three models. This means that at a given time step, the shaded region is defined 
by the maximum and minimum value of a given quantity seen among the three merger models. In the bottom panel, the yellow and blue bands denote the 
typical ranges of Xqq observed for the Galaxy and ULIRGs, respectively (as compiled by Tacconi et al. 2008). Prior to the burst, the inspiralling discs have 
Xqq values comparable to the Galactic mean. Upon the merger, increased velocity dispersions and gas temperatures contribute to lowering Xqq. In the 
post-merger stage, differences in H2 abundances, CO abundances, and time for the gas to re-virialise contribute toward a large dispersion in Xqq values. 



together. Consequently, we see a large range of X -factors in the 
GMCs throughout the galaxy. 

Because the emission from the merger is irregular, it is pos- 
sible that by choosing a different centre, the results from Figure 5 
would change. To test this, we recentered the image on the peak of 
the velocity-integrated intensity. Doing this provides no substantial 
change in the results of Figure 5. 



5 DISCUSSION 

5.1 Observational Consequences of the Model 

We have presented a model in which Xqq in GMCs is depen- 
dent on the physical conditions within the clouds. When the surface 
densities, kinetic temperatures and velocity dispersions within the 
GMCs resemble those of observed clouds in the Galaxy, the result- 
ing Xqq factor is comparable to the observed Galactic mean value. 
In starbursts, while the surface densities of clouds are higher, this 
is offset by both larger velocity dispersions in the GMCs as well 
as larger gas temperatures. The increased linewidths represent the 



turbulent velocity dispersion in the merger, as well as the stellar 
potential. The increased gas temperatures owe to efficient coupling 
with the dust at the high densities encountered in a merger. A fun- 
damental point of this study is that the physical conditions which 
cause Xqq to vary in starbursts are coupled. The same processes 
which drive the increased gas surface density also cause an increase 
in star formation rate which drives up the dust and consequently the 
gas temperatures. Similarly, in a merger-driven burst, the gas veloc- 
ity dispersion rises during the merger. 

While it is of utmost importance to parameterise Xqq in terms 
of observable properties, because the physical parameters which 
drive observed values of Xqq are coupled, this is a nontrivial task 
which is outside the scope of this work (though it will be inves- 
tigated in a forthcoming paper). Empirically, there is a tentative 
trend that Xqq decreases with increasing galaxy surface density 
(Tacconi et al. 2008). In the context of the models presented here, 
such a trend is plausible (see also Shetty et al. 2011b). One might 
expect that higher surface density systems typically arise in situ- 
ations when the velocity dispersion is high and the star formation 
rates, dust temperatures and gas temperatures are also high. 
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Figure 4. Maps of Xco f° r the fiducial disc galaxy and merger. The colour scale denotes logio (cm~ 2 /K-km s _1 ). Xco i s lower toward the centre of the 
disc galaxy due to higher temperatures and velocity dispersions in the clouds. The warm and high-cr gas is somewhat more spread out in the merger. We 
enforce a maximum Xco m the color bar of the disc galaxy of 21.3 to aid in clarity, though there are a few pixels with values as high as 21.8. 



While the models investigated here by no means comprise an 
exhaustive parameter-space study of galaxy masses, merger mass 
ratios, or merger orbits, we can investigate whether Xco can be 
parameterised by Eh2 in the simulations. To increase the dynamic 
range of surface densities in our models, we include one additional 
simulation of a high-redshift merger. The merger is the model sub- 
millimetre galaxy of Hayward et al. (201 1) during the coalescence 
when the peak merger-induced starburst is ~ 4500 IViQyr - . The 
model submillimetre galaxy has been shown to reproduce both the 
observed SED (Narayanan et al. 2010b), CO properties (Narayanan 
et al. 2009), overlap with 24 \im sources (Narayanan et al. 2010a), 
and number counts of observed SMGs (Hayward et al. 2010). Sim- 
ilarly, to increase the number of galaxies in our sample, we include 
many snapshots for the mergers (i.e. not just the snapshots at peak 
SFR). 

We plot the emission-weighted mean Xco of the GMCs in 
our model galaxies versus their emission-weighted mean surface 
densities in Figure 6. The models include the model mergers de- 
scribed in § 2, the model disc galaxy, and the newly introduced 
model submillimetre galaxy. The galaxies are binned by Eh2- The 
trend seen is what is expected: that Xco should decrease with gas 
surface density. In principle this owes to the fact that the highest 
surface density galaxies in our simulations are also forming stars 
at 10 2 — 10 3 MQyr -1 , and have relatively large velocity disper- 
sions. However, we emphasise strongly that Figure 6 is to be taken 
as a qualitative trend, rather than robust. A larger parameter-space 
survey of the simulations will be undertaken for a future study to 
identify the true mean and dispersions of Xco as a function of 

Eh2- 

The variation in Xco with environment may have implica- 
tions for observed Kennicutt-Schmidt star formation laws. An ex- 
ample of this was presented by Daddi et al. (2010) and Genzel et al. 
(2010) who applied a starburst Xco value to the inferred mergers 
in their observed sample of galaxies, and a Galactic Xco value 
to the discs. Doing so results in a bimodal star formation rate sur- 
face density-gas surface density relation. On the other hand, Os- 
triker & Shetty (2011) pointed out that if a simple E-dependent 



Xcooc E~ ' 5 is used above Eh2 = 100 Msun/pc 2 , a unimodal 
empirical star formation relation results (with Esfr oc E 2 , consis- 
tent with theoretical expectations for self-regulated star formation 
in this regime). 

Interpreting results for high-redshift galaxies in the context of 
our model is complex. Our model advocates for lower Xco values 
in high surface density environments. However, galaxies at z ~ 2 
which have surface densities comparable to local ULIRGs may in 
fact be discs (e.g. Daddi et al. 2005; Genzel et al. 2006; Forster 
Schreiber et al. 2009; Hopkins et al. 2010). How X C o of high-z 
discs should scale in this model is unclear at present. While their 
surface densities are comparable to local ULIRGs, their velocity 
dispersions may not show the same enhancement as seen in our 
model mergers. However, their gas temperatures may be compara- 
ble to their dust temperatures if the densities are high enough. Work 
by the Chicago group (R. Feldmann et al. in prep.) is underway 
to investigate this. Either way, the fact that our results tentatively 
suggest a relationship between Xco and surface density implies a 
continuum in Xco values, rather than a bimodality. Thus the rela- 
tionship between our model results and the interpretation of high-z 
observations by Daddi et al. (2010) and Genzel et al. (2010) will 
depend on the distribution of surface densities in their observed 
galaxies, among other issues. We note, however, that our work, like 
that of Teyssier et al. (2010), is consistent with the idea that the 
observed behaviour in the star formation rates and Xco can be ex- 
plained without the need to invoke a volumetric star formation law 
that is different in discs and mergers. The change in Xco we see 
in our simulations occurs because of changes in the physical con- 
ditions of GMCs associated with the merger, and not because the 
underlying star formation law is different. 

Finally, the concepts presented in this paper are testable in 
the near future with ALMA. Our models suggest that high spatial 
resolution observations of nearby ULIRGs will display both large 
velocity dispersions in the CO gas, and larger brightness tempera- 
tures than those seen in observations of Galactic GMCs on a com- 
parable scale. We see this when comparing the panels of Figure 1. 
Some observational evidence for this already exists. Interferomet- 
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Figure 5. Xco as a function of spatial distribution for the molecular gas in 
model disc (top) and merger (bottom). The Xco is derived from the maps 
(Figure 4), and is the emission-weighted mean Xco m bins of galactocen- 
tric distance. The bars around the points represent the range of Xco values 
within a given distance bin. The molecular gas at the centre of the disc has 
systematically lower Xco values than the outer disc. In the merger, Xco 
shows a wide-range of values throughout the galaxy. 



ric surveys of the central regions of nearby ULIRGs show veloc- 
ity dispersions of hundreds of km s~ , and brightness temperatures 
of tens of Kelvin (e.g. Scoville et al. 1997; Downes & Solomon 
1998). Similarly, unresolved observations of starbursts have shown 
gas and dust temperatures in the range of 30-50 K, in agreement 
with the models presented here (Yao et al. 2003; Narayanan et al. 
2005; Leech et al. 2010; Mao et al. 2010; Miihle et al. 2011). 



5.2 Relationship to Other Models 

The seminal work of Maloney & Black (1988) investigated Xco 
in galaxies via subresolution models of GMCs in a disc-like con- 
figuration. These authours found that Xco would vary from the 
Galactic value in cases of high kinetic temperature, high velocity 
dispersion or low metallicity. While not simultaneously modeling 
any of these effects, this model identified some of the most impor- 
tant driving factors in setting the observed X-factor in clouds. 

A number of other studies have also investigated Xco in mod- 
els of giant molecular clouds. Early studies implemented ID radia- 
tive transfer calculations in spherical models of GMCs (e.g. Kut- 
ner & Leung 1985; Wall 2007). With the increase of computational 
power, 3D numerical studies of GMCs in evolution have recently 
become feasible. Recently, Glover et al. (2010) and Glover & Mac 
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Figure 6. Emission-weighted mean Xco m GMCs for each model galaxy 
versus their emission-weighted mean surface density. The points represent 
the mean values for all of the GMCs within individual galaxies, and the 
galaxies are binned by surface density. In order to increase our sample size, 
many snapshots (i.e. not just the peak SFR for the mergers), as well as 
a model submillimetre galaxy, were included in this plot. Higher surface 
density environments tend to correspond with merger-driven starbursts in 
our models, and hence larger values of a and Tk- In these cases, Xco 
tends to be lower than the standard Galactic value. 



Low (2011) modeled H2 and CO formation/destruction in mag- 
netohydrodynamic models of GMCs. These models were elabo- 
rated upon by Shetty et al. (2011a,b) who utilised radiative trans- 
fer calculations in combination with these MHD models to pro- 
duce bona fide observables from the model clouds. These authours 
found that model GMC with mean densities, column densities, 
temperatures, and velocity dispersions comparable to the Milky 
Way's clouds (n ~ 10 2 - 10 3 cirT 3 , Xh 2 ~ 10 21 - 10 22 cm^ 2 , 
T ~ 10 — 20 K, a ~ 1 - 6 km s" 1 ) had average X C o fac- 
tors of order 2 — 4 x 10 20 cm~ 2 /K-km s _1 , and were insensitive 
to detailed temperature and velocity distributions. When manually 
increasing the velocity dispersion and/or temperature of the GMC, 
the resulting Xco values fell by a factor of ~ 5, comparable to 
both observed starbursts and the model mergers in this paper. While 
the simulations in Shetty et al. do not model the physical processes 
which may simultaneously cause Nh2, c, and T to vary, these mod- 
els do confirm that when increasing a or T and considering the ra- 
diative transfer through clouds, one will observe a depressed Xco, 
as is inferred in ULIRGs. 

In this sense, the models of Shetty et al. are complementary 
to those presented here. Shetty et al.'s models resolve much of the 
physics and chemistry within GMCs, though they have no infor- 
mation regarding the external environment from the host galaxy 
and how it may affect the cloud. Our simulations describe the am- 
bient environment surrounding the model GMCs, though at best 
they resolve the surfaces of the clouds (and require some amount of 
subresolution techniques). That both sets of models are converging 
upon the same result from different directions is encouraging. The 
next step forward will be to fully couple galaxy evolution simula- 
tions with high-resolution models of GMCs with a grid of model 
GMCs. These efforts are underway and will be presented in due 
course (R. Feldmann et al. in prep; Narayanan & Shetty in prep.). 
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6 SUMMARY AND CONCLUSIONS 

Utilising a combination of hydrodynamic simulations of disc 
galaxy evolution and galaxy mergers, dust and molecular line ra- 
diative transfer calculations, we investigated the dependence of the 
CO-H2 conversion factor on galactic environment. Our main results 
follow: 

(i) Provided that GMCs are gravitationally bound, disc galaxies 
in the local Universe have relatively little influence on the physical 
properties of GMCs within them (outside the central ~ kpc). The 
velocity dispersions are typically dominated by internal processes 
to the GMC, and the temperatures are roughly constant at ~ 10 K, 
set by a balance of molecular/atomic line cooling and cosmic ray 
and grain photoelectric effect heating. In this situation, when the 
surface densities of GMCs are comparable to those in the Galaxy, 
Xco will be similar to the Galactic value of Xco~ 2 — 4 x 10 2() 
cm _2 /K-km s _1 . 

(ii) In galaxy mergers, the GMC physical properties are strongly 
affected by the galaxy environment. The rise in surface density in 
GMCs during the merger is offset by an increase in the velocity dis- 
persion coupled to a rise in the kinetic temperature of the gas caused 
by efficient dust-gas thermal exchange at high densities. The com- 
bination of increased velocity dispersion and kinetic temperature 
increases the CO intensity, and lowers the observed Xco from the 
Galactic value by a typical factor of ~ 2 — 10. 

(iii) There is a slight trend with galactocentric radius such that 
GMCs toward the centres of disc galaxies will have a lower Xco 
than the disc-averaged value, owing to both increased velocity dis- 
persions in the clouds, as well as higher kinetic temperatures. 
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APPENDIX A: EFFECTS OF PARAMETER CHOICES 
AND ASSUMPTIONS 

In § 2, we outlined a number of parameter choices which could 
potentially influence the results in this paper. Here, we discuss the 
results in the context of these assumptions. 

Al Self-Consistency of The Temperature Calculations 

First, there is a discrepancy between the way the dust temperature is 
calculated in SUNRISE and in our temperature equilibrium model. 
In the former, the dust temperature is assumed to be in equilibrium 
with the radiation field, but we do not take into account any thermal 
exchange between the gas and dust. In the temperature equilibrium 
model, the dust grains are assumed to be able to exchange energy 
with the gas, but we hold the ambient radiation field fixed, rather 
than allowing it to change as the dust temperature does. Given the 
importance of the dust temperature in raising the gas temperature 
in this model, it is worth investigating any potential differences be- 
tween the two dust temperatures. 

In Figure Al, we plot the ratio of T dust from SUNRISE com- 
pared to T dust from the temperature equilibrium model as a func- 
tion of GMC density 9 for our fiducial merger. There is generally 
good agreement between the two, though some number of points 
at higher densities deviate strongly from unity. The gas which has 
poor agreement between the two dust temperature calculations al- 
most exclusively has all of its carbon in atomic form, rather than 
molecular. Beyond this, this gas tends to be towards the outskirts of 
the galaxy, in rather large cells in the adaptive mesh with relatively 
low masses (~ 10 Mq ). Because we enforce a rule that clouds 
must have a minimum surface density of 100 Mq pc~ 2 , these re- 
gions have extremely high densities, even if relatively low mass. In 
atomic gas of this density, the gas couples with the dust and can 
cause the dust temperature to change from that of the background 
radiation field. These outlying points have little effect on the final 
results, however, as they contain rather little mass. We denote the 
95th mass percentile by the blue shaded region in Figure Al. That 
is, the sum of the mass in the points outside of the blue region ac- 
counts for < 5% of the total molecular mass in the galaxy. As is 
shown, the differences in the two dust temperatures are small in this 
shaded region. 

A2 sunrise Input Parameters 

Similarly, a number of our assumptions in the SUNRISE modeling 
can have an effect on the derived dust temperature. We investigate 
those here. The important figure for these tests is Figure A2. Re- 
ferring to § 2.3, the time-averaged covering fraction of birthclouds 
around stellar clusters is a free parameter. While we chose a mod- 
est covering fraction, it is possible that a larger fraction may be 
reasonable. For example, if we choose a fraction /pdr = 1, the O 
and B stars would be blanketed by the ISM for their entire lives. 
This situation may exist in ULIRGs. To test for the effect in vary- 
ing /pdr , we plot the ratio of the dust temperatures derived from 
SUNRISE for a model with /pdr = 1 over the dust temperature 
from our fiducial model with /pdr = 0.3 against GMC density. 
The blue shaded region shows the 95th mass percentile, and grey 



9 This is not the actual mean density of the GMC, but the density account- 
ing for an enhancement by the turbulent compression of gas. This is the 
density that is used in the temperature equilibrium calculation. 
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Figure Al. Comparison between the SUNRISE dust temperatures and those 
calculated from the temperature equilibrium calculation as a function of 
density for model GMCs in our fiducial merger. See text for details. 



denotes the 99th mass percentile. The model with a larger clear- 
ing timescale for the birthclouds has cooler dust temperatures as 
less UV flux interacts with the ISM. However, the differences in 
dust temperature are generally within a factor of 50%, and much 
less than that (< 10%) when considering the bulk of the mass of 
the galaxy. We conclude that the PDR covering fraction is not an 
important driver in our model results. 

We can explore the effect of discarding the PDR birthcloud 
model, and assuming the cold ISM has a uniform volume filling 
fraction. In this case, the UV photons escape the star particles eas- 
ily, though optical depths for the photons in the ISM are large. For 
some number of the clouds outside of the nuclear region, the dust 
temperatures are thus colder than in our fiducial model (the low 
ratio points in Figure A2). However, the bulk of the gas mass is 
in a confined nuclear region in the galaxy which sees the intense 
UV radiation field. Because of this, this dust is heated well, and has 
comparable dust temperatures to our fiducial model. The dispersion 
in dust temperatures is again within 50%, and the 95th% percentile 
of mass shows relatively small discrepancies. 

It is possible that the AGN in the model merger contributes 
strongly to the dust temperature. To test this, we investigate a model 
where we have turned off the contribution of the AGN in determin- 
ing the dust temperature in the merger. While it is difficult to see 
by eye, the dust temperatures in the model with the AGN are hotter 
than the model with no AGN. By and large, however, the AGN is 
not powerful enough to have a significant effect on the overall tem- 
perature structure of the cold ISM as noted by the blue and grey 
shaded regions. 



A3 gadget-3 Input Parameters 

We now turn to possible parameters in the hydrodynamics simula- 
tions which may affect our results. As described in § 2.1, we make 
a number of parameter choices which may affect the star formation 
history of the model galaxies. Because the gas densities and dust 
and gas temperature can depend on these assumptions, it is worth 
exploring the robustness of our model results in the context of these 
choices. 
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Figure A2. Ratio of Zdust under various subresolution ISM specifications 
to that derived in our fiducial model versus H2 density in model GMCs. See 
text for details. 



There are two principle parameter choices which govern the 
physical state of the ISM in our hydrodynamic modes: the star for- 
mation "law", and the equation of state. As discussed in § 2.1, we 
adopt a star formation law such that the star formation time scale 
is assumed to be proportional to the local dynamical timescale, and 
whose rate matches the normalisation of the locally observed Ken- 
nicutt (1998a) relation. 

In the absence of a complete theory of star formation, a num- 
ber of possible choices exist regarding the implementation of a star 
formation recipe on subresolution scales. One can imagine a similar 
Kennicutt-Schmidt solution, though with an index of unity as ap- 
pears to be observed on resolved scales in nearby galaxies (Bigiel 
et al. 2008), or steeper index as tentative observational (Bouche 
et al. 2007) and theoretical (Feldmann et al. 2011) evidence mo- 
tivates at high-redshift. Beyond this, more sophisticated physical 
models may provide reasonable prescriptions for star formation in 
galaxy evolution models (e.g. Tan 2000; Krumholz et al. 2009b; 
Ostriker et al. 2010). While it is outside the scope of this study 
to perform a detailed study of various star formation recipes in 
GADGET-3, we perform some simple tests to investigate the role 
of our adopted star formation law in driving the simulated Xco 
factors. 

In the Table, we show the luminosity-weighted mean Xco 
values for our fiducial merger during the peak of their starburst for 
a variety of input parameters in our GADGET-3 simulations. The 
Kennicutt-Schmidt index=2 case shows a comparable Xco value 
as the fiducial merger. While the burst is moderately diminished 
(owing to rapid consumption of the gas during early phases of the 
merger's evolution, though still quite large at ~ 250 MQyr -1 as 
opposed to ~ 340 MQyr -1 in the fiducial merger), the large stellar 
mass upon coalescence maintains a large velocity dispersion in the 
gas. This drives a low mean Xco • 

The model with the largest impact on our results is the model 
with Kennicutt-Schmidt index=l; This model merger has a mean 
Xco comparable to the model disc. Models with a KS index of 1 
do not go through a burst upon merging (T. Cox, private communi- 
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cation). With a KS index of 1, to first order, the total star formation 
rate is proportional to the total gas mass. Because we don't include 
any gas replenishment from the intergalactic medium, the gas mass 
only decreases with time, as does the SFR in this simulation. The 
low SFR upon merging leads to low gas/dust temperatures, and in- 
creased Xco- We note that this situation is unlikely to describe 
real mergers as, observationally, galaxy mergers exhibit the high- 
est star formation rates in the local Universe. Both observational 
(Bigiel et al. 2008) and theoretical (Krumholz et al. 2009b; Os- 
triker & Shetty 2011) evidence suggest that a KS index > 1 may 
describe high-surface density systems. 

We utilise an equation of state for the ISM which incorpo- 
rates a subresolution prescription for capturing the effect of super- 
novae heating of the ISM (Springel & Hernquist 2003). The nom- 
inal Springel & Hernquist (2003) and Springel et al. (2005b) EOS 
is given by: 

fcff = (7 - l)(PhM h + PcUo) (Al) 

where 7 = 5/3 as the adiabatic index of the gas, ph,c is the density 
of the hot and cold phase, and «h,c is the specific thermal energy of 
the two phases. For a given IMF, Springel & Hernquist (2003) show 
that the EOS is completely defined by the star formation time scale, 
the normalisation of the cloud evaporation rate, and a supernovae 
"temperature" which defines the heating rate from supernovae of 
a given IMF. The full "effective" EOS has the property in which 
pressure rises with density faster than an isothermal gas, as can be 
seen in Figure 4 of Springel et al. (2005b). 

Our fiducial model utilises a softer EOS than the full model. In 
particular, we interpolate between the full "stiff" model (where we 
assign a parameter <jeos = 1) and isothermal model (<?eos = 0) 
and employ qeos = 0.25. In the full "stiff" EOS as in simulations 
of discs scaled for the local Universe, the ISM can become so pres- 
surised as to appear smooth with relatively few clumps (Springel 
et al. 2005b). To test our assumption of a softer EOS, we have run 
a test simulation with g EO s=l. 

In the Table, we show the mean Xco for our fiducial merger, 
though with g EO s=l (denoted "Stiff EOS"). We see a larger value 
than the fiducial model. Isolating the root cause is nontrivial. By ef- 
fectively increasing the effect of supernovae feedback, we increase 
dust heating (by reducing the dumpiness of the gas), though we 
also reduce the magnitude of the burst (due to a retardation of gas 
fragmentation). These effects serve to somewhat offset one another 
with respect to the gas temperature. 

We test whether our assumption of Go = 1 outside of clouds 
plays a strong role in our model results. Following Ostriker et al. 
(2010), we have run a model in which we scale the interstellar ra- 
diation field by the value of the local SFR compared to that in the 
solar neighbourhood, and show the mean Xco value in the Table 
(denoted by "OML G-Scaling"). Because the clouds are strongly 
shielded in this model, scaling G makes little difference. 

Finally, we consider the spatial resolution of our model. Our 
current resolution has cell sizes (~ 70 pc) which are of order 
the SPH smoothing length. Further increasing the spatial resolu- 
tion does not provide new physical information. More seriously, 
increasing the spatial resolution of the SPH simulations would 
run into scenarios of unphysical descriptions of the ISM with the 
Springel & Hernquist (2003) multiphase model. We can see, how- 
ever, what direction the results may go if it were possible to in- 
crease our resolution. To investigate this, we have run our fiducial 
merger snapshot with one less level of refinement in the adaptive 
mesh, giving a minimum cell size of ~ 140pc, and show the mean 
Xco in the Table (denoted by "Half-Res"). The mean Xco in the 



Table Al. 



Model 


Mean Xqq 




(cm^/K-kms- 1 ) 


Fiducial 


6.31 X 10 19 


Stiff EOS 


1.08 x 10 20 


KS = 2 


5.43 x 10 19 


KS = 1 


2.73 x 10 20 


OML G-scaling 


7.10 x 10 19 


Half-Res 


3.02 x 10 19 



low resolution model is 60% the value of the fiducial model. It 
is conceivable, then, that further increasing the spatial resolution 
would increase Xco in the merger model. On the other hand, re- 
ducing the cell size would result in more resolved GMCs. Because 
it is the resolved GMCs that drive the mean Xco for the merger 
model, it is also possible that increased spatial resolution would 
result in little change to the results presented here. 



